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ABSTRACT 

Wc derive a model-independent expression for the minimum line-of-sight extent of the hot plasma in 
a cluster of galaxies. The only inputs are the 1-5 keV X-ray surface brightness and the Comptonization 
from Sunyaev-Zel'dovich (SZ) data. No a priori assumptions regarding equilibrium or geometry are 
required. The method applies when the X-ray emitting material has temperatures anywhere between 
0.3 keV and 20 keV and metallicities between and twice solar — conditions fulfilled by nearly all 
intracluster plasma. Using this method, joint APEX-SZ and Chandra X-ray Observatory data on the 
Bullet Cluster yield a lower limit of 400 ± 56 kpc on the half-pressure depth of the main component, 
limiting it to being at least spherical, if not cigar-shaped primarily along the line of sight. 
Subject headings: Galaxies: clusters: general — cosmic microwave background — X-rays: galaxies: 
clusters — radiation mechanisms: thermal 



1. INTRODUCTION 

The hot intracluster medium in clusters of galaxies 
emits X-rays via thermal Bremsstrahlung and line emis- 
sion; the very same atmosphere blueshifts cosmic mi- 
crowave background photons via the Sunyaev-Zel'dovich 
(SZ) effect (for a review see Birkinshaw 1999; Carlstrom 
et al. 2002, and references therein). The X-ray emissivity 
scales as electron density squared, but the SZ spectrum 
modification depends on the integral of pressure along 
the line of sight (and hence only linearly on the electron 
density). As a result, if the cosmology is known, it is pos- 
sible to recover structural information about the cluster 
atmosphere (Myers et al. 1997; Hughes & Birkinshaw 
1998; Grego et al. 2000). Such "deprojection," however, 
has always required drastic assumptions regarding the 
3D structure of the ICM; in general, imposing hydro- 
static equilibrium is often required (Zaroubi et al. 2001). 
Most importantly, spherical or triaxial geometry have al- 
ways been necessary assumptions as well (Pointecouteau 
et al. 2002; Lee & Suto 2004; De Filippis et al. 2005; 
Ameglio et al. 2007, 2009; Nord et al. 2009; Basu et al. 
2010; Allison et al. 2011). So far, forward convolution of 
simple, smooth, parmctcrizcd ICM thermodynamic vari- 
ables has been the only way to tie the X-ray and SZ 
observables together to determine cluster geometry. 

However, a large fraction of clusters are not relaxed 
and are poorly described by hydrostatic or spherically 
symmetric models (e.g. Eckert et al. 2011). The most 
spectacular and interesting of these systems are violent 
mergers such as the Bullet Cluster (Clowe et al. 2006), 
Abell 520 (Mahdavi et al. 2007), or MACSJ0025 (Bradac 
et al. 2008). These systems clearly resist any attempt 
at description using symmetric equilibrium models. At 
the same time, it would be highly valuable to recover 
some information regarding the structure of such objects 
along the line of sight. Such information could be used to 
constrain the parameter space of supercomputer N-body 
models that seek to replicate and thus better understand 
the collisions (e.g. Mastropietro & Burkert 2008; Randall 
et al. 2008; Forero-Romero et al. 2010). 



In this Letter, wc show that the Cauchy-Schwarz inte- 
gral inequality can yield useful information on the depth 
of a chister of galaxic;s without assumptions regarding 
the geometry or thermodynamic state of the ICM. In §2 
we derive the inequality, and in §3 we interpret it. In §4 
we apply it to the Bullet Cluster as a test case, and in 
§5 we conclude our discussion. 

2. THEORETICAL FRAMEWORK 

In a cluster of galaxies, the outgoing surface brightness 
of a column of hot, single-phase X-ray emitting plasma 
is 

S = / dE nenHAEiT,Z)dz (1) 

Jei J-oo 

where the emission is in a rest-frame bandpass between 
photon energies Ei and i?2- The electron and hydrogen 
densities are rig and uh, T is the plasma temperature, 
Z is the plasma metallicity in solar units, and A^; is 
the cooling function incorporating the bremsstrahlung 
continuum as well as line emission. For convenience we 
can recast this formula as 

/oo 
nlABp{T, Z)dz (2) 
-oo 

where Abp (T, Z) is the cooling function integrated over 
the X-ray bandpass and multiplied by the electron-to- 
proton ratio, {ne/nn) of the plasma, which varies less 
than 1% for typical ICM metallicities. 

In the non-relativistic limit, the fractional Sunyaev- 
Zel'dovich decrement is proportional to the Compton y 
parameter: 

1 /"^ 

y = ^ / TleOTkTdz (3) 

where gt is the Thomson cross section, rrie is the elec- 
tron mass, c is the speed of light, and k is Boltzmann's 
constant. We now also define the half-pressure depth zr 
such that 

2 f^" 

y = / ^eCTTkTdz (4) 
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where the origin of the 76 axis is arbitrary so long as 
the above equation is fulfilled. In other words, half the 
integrated pressure of the cluster along the line of sight 
is confined to a region of size zh- 

The X-ray surface brightness and Compton y param- 
eter involve integrals of the thermodynamic variables 
along the line of sight. They therefore cannot be inverted 
without assumptions regarding the geometry of the clus- 
ter. However, consider the Cauchy-Schwarz inequality: 

y" f{x)g{x)dx < J fixfdx J g{xfdx (5) 

Here we show that this integral is useful in recasting 
the depth measurement problem in model-independent 
terms. Specifically, consider the choice / = n^T, g = 1; 
then, integrating only along the half-pressure depth, we 
have 

dz > H° , , 



r 

JO 



(6) 



Furthermore, since the square of any function is more 
compact that the function itself, equations (3) and (4) 
together imply that 



iT^dz > 2 



tT^dz 



(7) 



The above inequality follows from the fact that the vari- 
ance of any square-integrable function /(x) is equal to or 
larger than the variance of f{x)'^. Therefore it is always 
possible to choose a z-axis origin such that [0, zh] con- 
tains half the integrated pressure but more than half the 
integral of the pressure squared. We therefore obtain 



" dz > 2l&^ 



(8) 



Combining equation (8) with the Comptonization pa- 
rameter (equation 3), and multiplying by S/E = 1, we 
can write 

1 J — OO *2 

This constraint on cluster depth can be rewritten as 

ZH>£xC (10) 
where we define i as the fiducial depth of the cluster 
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^0.56 



lO-*''^ erg s-i Mpc" 



(11) 



Mpc(12) 



and where we define the dimensionless, bandpass- 
dependent thermodynamic uncertainty factor 



TTZeC / n\h.'QY>{T, Z)dz 

aTjnlk^T^dz 



(13) 



If we were dealing with a medium of constant tempera- 
ture and metallicity, we could now pull out the functions 
of T and Z, cancel the emission measure integrals, and be 



done; hut we do not make this assumption and continue 
with the general case. 

So far we have shown that given a measurement of the 
Compton y parameter and the X-ray surface brightness, 
the half-pressure depth of the cluster along the line of 
sight is at least £ (equation 12) times a correction factor 
C (equation 13). 

We now attempt to derive useful bounds on C. Sup- 
pose that, for any given X-ray bandpass, ( > Clow for all 
plausible values of T and Z. In that case without any 
loss of generality we can write 



Zh>£x Clow 



(14) 



To find Clow 5 imagine three arbitrary one dimen- 
sional positive functions ne{z), p{z) and q{z). If 
p/q > Clow everywhere, then necessarily p > Ciow9i 
and so n^p > Ciow'T^e? everywhere. It follows that 
J^ne{z)'^p{z) > Clow /a '^e(2)^a(-2)- In other words, the 
ratio of two one-dimensional integrals is always greater 
than the smallest value taken on by the ratio of the in- 
tegrands, as long as the integrands are positive and the 
limits of integration are the same. 

Choosing p = mecABp(T, Z) and q = ark'^T'^, it fol- 
lows that C > Clow if 



Ci 



low 



minimum 



of 



mecABp(T,Z) 



(15) 



for all Z and T of interest. Our task then becomes one 
of minimizing the quantity 

mecABp(T,Z) 
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(16) 



3. INTERPRETATION 



The lower limit on the half-pressure depth of a clus- 
ter of galaxies becomes £ x Ciow, where Ciow is the min- 
imum of the bandpass-dependent quantity 7 over all 
temperatures and metallicities contributing to the SZ 
and X-ray fiuxes. Since the cooling function for ther- 
mal Bremsstrahlung A(r, Z) ~ for T » 10 keV, it 
follows that 7 ~ T~^/2 for large T, and the Zh con- 
straint requires the assumption of a high-temperature 
cutoff in order to be relevant (i.e., to avoid the trivial re- 
sult Zh > 0). We call this high-temperature cutoff Tmax- 
We leave Tmax flexible on a cluster-by-cluster basis. Its 
value can be constrained via modeling of sufficiently high 
signal-to-noise X-ray spectra, or else set to a sufficiently 
high value (e.g. no cluster is believed to possess thermal 
gas beyond 20 keV). 

A low temperature cutoff, Tj^in, is also needed, be- 
cause Clow — whenever Abp(T, Z) 0, i.e., when 
there is a column of material emitting outside the X- 
ray bandpass. In a bandpass centered on photon en- 
ergy E, all emission with T > E/3 contributes to the 
X-ray flux, but emission with T < E/3 contributes neg- 
ligibly. For most common X-ray bandpasses, this trans- 
lates to T < 0.3 keV. Searches for very cool, T < 0.3 
keV material — typically referred to as the warm-hot in- 
tergalactic medium or WHIM — have thus far been in- 
conclusive. While the WHIM must exist, it likely has 
densities < 5 x 10^''' cm~^ (Nicastro et al. 2005), and 
thus would contribute negligibly to the overall SZ and 
X-ray flux in the direction of a massive cluster, where the 
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typical gas densities are 10^* cm^'^ or higher. We there- 
fore assume that in any source with both a significant 
SZ decrement and a measured Chandra/XMM-Newton 
X-ray flux, the contribution of material at T < 0.3 keV 
is negligible. For example, consider a medium at T = 0.3 
keV, He = 5 X 10~^ cm~^, emitting over a length of 10 
Mpc along the line of sight. Such a medium would make 
an additive contribution of 6 x 10^^ to the central Comp- 
ton y parameter, a less than one percent effect for massive 
clusters even in the maximal WHIM case. We therefore 
neglect the WHIM for our purposes and set Tmin = 0.3 
keV. 

As the bandpass to focus on in the present Letter, we 
select the restframe 1-5 keV photon energy range. This 
energy range has the advantage of overlapping the most 
sensitive regions of current X-ray telescopes, and of be- 
ing detectable between redshifts z = Q and z = 1 with 
relatively low astrophysical background contamination. 

With the bandpass chosen, we numerically minimize 7 
as a function of temperature and metallicity, assuming a 
MEKAL (Mewe et al. 1985) plasma with allowed metal- 
licities from to twice solar, and allowed temperatures 
ranging from Tmin to Tmax- We find that the minimum 
value of 7 is sensitive only to T^ax- For all Tnjjn > 0.3 
keV, Tmax > 3.5 keV we find that the following empirical 
formula fits the exact value of Ciow with better than 1% 
accuracy: 



Clow = 1-883 + 



54.12 



19.05 



(17) 



I-t!- Clow has no dependence on the coolest material so 
long as its temperature is above 0.3 keV. We note that 
beyond T — 2Q keV the relativistic correction to equa- 
tion (3), which we neglect here, becomes important. The 
above formula is therefore only valid between T^ax = 3.5 
and 20 keV. 

4. APPLICATION TO THE BULLET CLUSTER 

To demonstrate this technique, we jointly analyze X- 
ray and Sunyaev-Zel'dovich observations of the Bullet 
Cluster. This cluster is a particularly apt candidate for 
this type of analysis, due to the clear incorrectness of 
typical spherically symmetric models in describing its 
morphology. The Bullet Cluster consists of a "main" 
diffuse component plus a "bullet" to the west. Us- 
ing 150GHz APEX-SZ, Halverson et al. (2009) infer a 
Comptonization for the main component of the cluster 
of j/o = 3.4±0.3x 10^*. This Compton y measurement is 
not model-independent — it relies to some extent on the 
precise 2D shape of the cluster in the plane of the sky, 
reflecting the limited angular resolution of APEX-SZ ex- 
periment. Further uncertainties arise from the fact that 
some leakage from the Bullet 1.5' to the west into the cen- 
tral APEX-SZ beam is likely. It is difficult to estimate 
the bias in y introduced here, first because the Bullet is 
both denser and cooler than the main component of the 
cluster, and second because the Compton y model adopts 
the X-ray profile as a strong prior. A further issue is that 
the azimuthally averaged profile has better statistics at 
radii which include the Bullet (and therefore the fit may 
be particularly driven by those regions). These short- 
comings highlight the need for higher angular resolution 
SZ experiments if bias-free 3D measurements of cluster 



structure are to be achieved. Such observations would 
allow for more complicated 2D Compton y distributions 
and thus overcome much of the systematic problems of 
the smooth symmetric models. 

To measure the rest-frame 1-5 keV surface brightness, 
we reanalyze an archived lOOks Chandra X-ray Observa- 
tory data set (ObsID 5356). We extract a region of radius 
182 kpc (0.691' at Hq = 71 kms'^ Mpc-\ = 0.23) 
centered on the main portion (i.e., the true X-ray center 
of the difi^use component of the cluster). The aperture 
was selected to match the beam size of the APEX-SZ ex- 
periment. We refer to Mahdavi et al. (2007) for all X-ray 
data reduction procedures, including particle background 
subtraction, point source masking, and other details. 

Using a redshift z = 0.296, we obtain a best-fit temper- 
ature of 11.5 ± 0.6 keV, best-fit metallicity of 0.23 ± 0.05 
solar, and best-fit absorbing column of 7.1 ± 0.5 x 10^° 
cm~^. The X-ray flux in the 0.77-3.9 keV band (cor- 
responding to the 1-5 keV rest frame band) is 1.89 ± 
0.01 X 10~^^ ergs cm~^ s~^. Multiplying by Att times the 
square of the luminosity distance at the above cosmol- 
ogy (1517 Mpc), and dividing by the projected area of 
0.104 Mpc^, we obtain a restframe surface brightness of 
5.03 ± 0.03 X 10*5 ergs s'^Mpc-^. 

Via equation (12) the observations therefore imply an 
upper limit 



ZH > [0.56 X (3.2)2 ^ (5.0)-!] (low 



(18) 



To constrain Ciow) we need an estimate of the highest tem- 
perature gas along the line of sight. To be conservative, 
we assume that temperatures up to 20 keV contribute 

signiflcantly to the emission measure (this is consistent 
with the measured X-ray spectrum). Plugging 20 keV 
into equations 17 and 18 yields a lower limit 



zh > 400 ± 56 kpc 



(19) 



where the dominant source of error is the statistical un- 
certainty on the Compton y parameter. The (unknown) 
systematic bias in y is not included. 

Now we attempt to calculate the axial ratio of the main 
(eastern) component of the cluster, comparing the line- 
of-sight half-pressure depth zh to the plane-of-sky half- 
pressure width, which we call wh- Halverson et al. (2009) 
flnd that 

AT/T <x{l + 6»2/6'^)i-3'3/2 (20) 

where /3 = 1.2 ± 0.13 and Oc = 142" ± 18". Numerically 
integrating this profile, we find that half the flux is con- 
tained within a circular aperture of radius 0.679c, which 
translates into a half-width width of ss 2 x 0.67 x 325 
kpc ss 435 kpc. Note that this is likely more an upper 
limit than an accurate measurement, given the mixing 
of the bullet and main component signals due to poor 
resolution of the APEX-SZ instrument. However, the di- 
rection of this inequality works in our favor, because a 
lower limit on zjf divided by an upper limit on Wh still 
yields a valid lower limit on zh/wh- 



Zh 
Wh 



> 0.92 



(21) 



Thus, the ratio of the half-pressure diameter of the 
Bullet Cluster to its half-pressure width in the sky is 
constrained to be greater than approximately one. 
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Fig. 1. The dimcnsionlcss thermodynamic uncertainty factor 
(^low '^^ a function of maximum plasma temperature Tmax, calcu- 
lated for the restframe 1-5 keV bandpass. The curve is given by 
equation 17. 



5. CONCLUSION 

Wc derive the first completely model- and geometry- 
independent expression limiting the half-pressure depth 
of a galaxy cluster along the line of sight. Regard- 
less of geometry, clumping, or thermodynamic state of 
the intracluster medium, combining X-ray and Sunyaev- 
Zel'dovich observations of a single cluster can yield a 
lower limit on this depth. This lower limit is 



m^c'y 



(22) 



where Ciow is a dimensionless number which depends on 
the specific bandpass used and the maximum ICM tem- 
perature that contributes significantly to the SZ effect. 



Equation (17) gives an empirical fitting formula for C,\o^ 
in the 1-5 keV restframe bandpass (also shown in Figure 
l). 

This lower limit is a function of largely model- 
independent observables: y (the Compton y parameter), 
and S (the X-ray surface brightness). No deprojection, 
and no modeling beyond standard reduction of the SZ 
decrement and X-ray spectrum is required. The method 
neglects relativistic SZ corrections, and assumes that the 
bulk of the pressure and cooling is done by gas with tem- 
peratures between 0.3 keV and 20 keV, and metallicities 
between and 2 times solar. Otherwise, the gas can have 
an arbitrary distribution along the line of sight. The 
constraint can be made more stringent by relaxing the 
upper limit of gas temperatures allowed; one way to do 
this in the future is via measurements of the differential 
emission measure dn^/dT. This requires some more so- 
phisticated spectral modeling, but still no assumptions 
regarding geometry or equilibrium. 

These results are useful in a number of scenarios. For 
merging clusters of galaxies such as the Bullet Cluster, 
they can help constrain the geometry of the merger. Us- 
ing this method, our analysis of joint APEX-SZ and 
Chandra data for the Bullet Cluster imply a minimum 
line-of-sight half-pressure depth of > 400 kpc, implying 
a minimum line-of-sight to plane of the sky axial ratio of 
« 1. 

Thus the Bullet Cluster is constrained to be at least 
spherical, if not cigar-shaped primarily along the line of 
sight. For larger samples of X-ray emitting clusters, this 
technique can help quantify possible selection biases to- 
wards line-of-sight elongations, as well as possibly aid in 
determining the geometry of shock and cold fronts. 
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